Revision of model parameters for K-type charge transfer salts: an ab initio study. 
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Intense experimental and theoretical studies have demonstrated that the anisotropic triangular 
lattice as realized in the k-(BEDT-TTF)2X family of organic charge transfer (CT) salts yields 
a complex phase diagram with magnetic, superconducting, Mott insulating and even spin liquid 
phases. With extensive density functional theory (DFT) calculations we refresh the link between 
manybody theory and experiment by determining hopping parameters of the underlying Hubbard 
model. This leads us to revise the widely used semiempirical parameters in the direction of less 
frustrated, more anisotropic triangular lattices. The implications of these results on the systems' 
description are discussed. 

PACS numbers: 74.70.Kn,71.10.Fd,71.15.Mb,71.20.Rv 
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A strong research trend of the new millennium has 
been the desire to understand complex manybody phe- 
nomena like superconductivity and magnetism by realis- 
tic modelling, i.e. to employ precise first principles calcu- 
lations to feed the intricate details of real materials into 
the parameter sets of model Hamiltonians that are then 
solved with increasingly powerful manybody techniques. 
The K-(BEDT-TTF)2X [l| organic charge transfer sahs 
are a perfect example for a class of materials with such 
fascinating properties that they drive progress in exper- 
imental and manybody methods alike. Experimentally, 
the phase diagram shows Mott insulating, superconduct- 
ing, magnetic and spin liquid phases [H, d, |j, Q . Theo- 
retically, the underlying anisotropic triangular lattice is 
a great challenge due to effects of frustration and the 
intense efforts to get a grip on the problem include stud- 
ies with path integral renormalization group (PIRG) d| , 
exact diagonalization 7^, variational Monte Carlo 8|, 
cluster dynamical mean field theory and dual 

Fermions [Tlj to cite a few. In this rapidly expanding 
field of research, electronic structure calculations play 
the decisive role of mediating between the complex un- 
derlying structure and phenomenology of organic charge 
transfer salts and the models used for understanding the 
physics jT2|, and in this work, we will provide the per- 
spective of precise, state of the art electronic structure 
calculations. 

Previously, K-type CT salts have been investigated by 
semi-empiricial and first principles electronic structure 
calculations. The most commonly used t, t' , U parame- 
ter sets derive from extended Hiickel molecular orbitals 



calculations [13|, llj| performed on different constellations 
of BEDT-TTF dimers. 

The main result reported here is that our first prin- 
ciples study shows all four considered K-type CT salts 
to be less frustrated than previously assumed based on 
semiempirical theory. Most importantly, the often cited 
value of t' /t = 1.06 [l3| for the spin liquid material k- 
(ET)2Cu2(CN)3 should be replaced by the significantly 



smaller value t'/t = 0.83 ± 0.08. This has fundamen- 
tal implications on the systems' model description as we 
shall see below. 

In this Letter, we emplo y t he Car-Parrinello [l3| 
projector-augmented wave [iG] molecular dynamics 
(CPMD) method for relaxing the only partially known 
structures and per form structure optimizations at con- 
stant pressure [l7| in order to prepare high pressure 
structures. We calculated the electronic structure us- 
ing two full potential all-electron codes, the linearized 
augmented plane wave (LAPW) method implemented in 
WIEN2K[i|] and the full potential local-orbital (FPLO) 
method [l9|. 

We base our study on the following structures: For k- 
(BEDT-TTF)2Cu2(CN)3, the crystal structure ^ was 
published without hydrogen atom positions; we add them 
and perform a CPMD relaxation within the general- 
ized gradient approximation (GGA) in which we keep 
lattice parameters and heavy atom (S, Cu) positions 



fixed [21|. With this procedure we make sure that the 



important atomic positions in the crystal structure re- 
main as given by the diffraction measurements. We 
remove the inversion center in the middle of a cyano 
group by lowering the symmetry from P2i/c to Pc. No 
high pressure structures of k-(BEDT-TTF)2Cu2(CN)3 
are published, and we therefore perform a full relax- 
ation of the structure and lattice parameters at a pres- 
sure P — 0.75 GPa [l^l; we recover the Pc symmetry 
by further optimizing the high pressure structure with 
symmetry constraints. For k-(BEDT-TTF)2Cu(SCN)2, 
structures at ambient pressure and at P = 0.75 GPa are 
known . Ambient pressure structures are available for 
K-(BEDT-TTF)2Cu[N(CN)2]Br [H, ^ and k-(BEDT- 
TTF)2Cu[N(CN)2]Cl Q. We will refer to the four ma- 
terials as K-CN, K-SCN, K-Br and k-CI for brevity. At 
ambient pressure and low temperatures k-SCN and K-Br 
superconduct while k-CI orders antiferromagnetically. 

As an example for the structure of all K-type CT salts 
considered here, we show in Fig.[l]the structure of k-CN. 



FIG. 1; (Color online) Structure of k-(BEDT- 
TTF)2Cu2(CN)3 as typical example of the K-type CT 
salts, (a) Side view with ET layers separated by insulating 
Cu2(CN)3 anion layers, (b) ET lattice, viewed in the be 
plane. The unit cell containing four ET molecules is shown, 
(c) ET lattice shown in the same projection as (b) where 
the ET molecules have been replaced by single spheres. ET 
dimers are emphasized by heavier bonds, (d) Lattice of ET 
dimers, replaced by single spheres. 
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FIG. 2: (Color online) DFT bandstructures of (a) k-(BEDT- 
TTF)2Cu2(CN)3 at ambient pressure and at P = 0.75 GPa, 
and of (b) k-(BEDT-TTF)2Cu[N(CN)2]C1 shown with sym- 
bols. Lines represent the tight binding fits to the bands close 
to Ef deriving from ET molecules. 



The strongly anisotropic nature of the materials is ap- 
parent from the alternation of charge donating ET layers 
and acceptor anion layers (Fig, [l] (a)), Cu2(CN)3 in this 
case. Experimental observations [iSj as well as our calcu- 
lations show that Cu occurs in Cu"''^ oxidation state here 
and has a filled 3d shell, making the anion layers insulat- 
ing. The K-type arrangement of ET molecules (Fig.[T](b)) 
exhibits a strong dimerization of the molecules, four of 
which are contained in a unit cell. They form the lattice 
shown in simplified form in (Fig. [1] (c)). In the charge 
transfer process, each molecule donates a charge of — 0.5e 
and is left with half a hole. Considering pairs of dimers 
as the fundamental unit leads to the triangular lattice of 
Fig. [1] (d) with one hole per dimer. 

In Fig. [2] we show three examples of low energy band- 
structures for K-CN at ambient and elevated pressure 
(Fig. El (a)) and for k-C\ (Fig. [2] (b)). The overall 
features of the bandstructures are very similar to the 
known semiempirical band structures, with two anti- 
bonding bands of the ET molecules crossing the Fermi 
level, and the corresponding bonding bands nearly mirror 
symmetric below. Due to the double number of eight ET 
molecules per unit cell in the case of k-C1, Fig. [2] (b) ac- 
tually has four bands crossing the Fermi level, which are 
pairwise degenerate due to the high symmetry (Pnma) 



of the crystal structure. In contrast to previous calcula- 
tions, there are four bands with mainly Cu 3d character 
close to the ET bands; slightly below the bonding ET 
bands in the case of k-CN and slightly below the anti- 
bonding ET bands in the case of k-C1. The bandstruc- 
tures of Fig. [5] (a) were performed with LAPW using the 
GGA functional 27|. We checked the reliability and re- 
producibility of these results by repeating the calculation 
with FPLO 29], and we found that the band structures 
and all derived quantities match very well. As k-C1 has 
roughly twice as many atoms in the unit cell compared 
to K-CN, we used the faster FPLO method, within GGA, 
for Fig. [21(b). The low energy bands show almost no dis- 
persion in the direction perpendicular to the anion layers 
which is the a direction (M-Y in the Brillouin zone) for 
K-CN and the b direction (U-Z in the Brillouin zone) for 
K-Cl, confirming the quasi-two-dimensional character of 
the K-type CT salts. 

We corroborated that the calculated bandstructures 
agree with experimental evidence by calculating the 
Fermi surface (FS) for k-CN at ambient pressure and 
at P = 0.75 GPa. The FS consists of two different sheets 
arising from the two bands crossing the Fermi energy, 
a corrugated cylinder around Z and a quasi one dimen- 
sional sheet parallel to M — Y. Both sheets touch along 
the line M — Z and give rise to a combined elliptical orbit 
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FIG. 3: (Color online) Fermi surface cuts for nonmagnetic 
K-(BEDT-TTF)2Cu2(CN)3 at ambient pressure (black) and 
at P = 0.75 GPa (light). Thin lines show the respective first 
Brillouin zone. The combined elliptical orbit as detected by 
AMRO measurements (see text) is indicated by broken lines. 

(Fig. [31) which encloses 99-100% of the area of the first 
Brillouin zone. This is in good agreement with previous 
calculations [13] and with angle-dependent niagnetore- 
sistance oscillations (AMRO) measurements f30'|. Our 
calculated FS shows the same behavior with pressure as 
observed in experiment. The area enclosed by the com- 
bined elliptical orbit increases by 6% from 37.4 nm~^ 
at ambient pressure to 39.7 nm~^ at P = 7.5 kbar (see 
Fig. which is consistent with the increase of 6% be- 
tween 2.1 and 7 kbar obtained from AMRO measure- 
ments in Ref. f30| . 

In the next step, we use the DFT band structures to 
extract the parameters of the Hubbard Hamiltonian 

H= J2 t(4s-+H.c.)+^i'(cLs,+H.c.) 

<ij>,cr [ij],cr 
i 

where < ij > and [ij] indicate sums over nearest and next 
nearest neighbors, respectively. We do this by fitting the 
bands to a tight binding model. Note that also nonlocal 
correlations V, V are thought to play a role in organic 
CT salts but their determination is beyond the scope 
of our present investigation. We consider each molecule 
position as shown in Fig. [1] (c) as a site; then the four 
(eight) sites per unit cell contribute the four (eight) ET 
bands at the Fermi level in the case of k-CN and k-SCN 
(K-Br and k-CV) This yields the hopping integrals ti to 
ti where the index indicates increasing distance; ti is the 
intradimer hopping integral. The black lines in Fig. [2] 
show as an example the tight-binding fit obtained from 
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FIG. 4: (Color online) Overview of hopping parameters ob- 
tained from GGA (filled symbols) and LDA (open sym- 
bols) calculations for k-(BEDT-TTF)2Cu2(CN)3, k-(BEDT- 
TTF)2Cu(SCN)2 and «:-(BEDT-TTF)2Cu[N(CN)2]C1. The 
spread of t and t' parameters depending on basis set, func- 
tional and number of fitted bands yields t' /t ratios with a 
margin of error indicated by lines and shaded regions. Ar- 
rows indicate the change of t' /t ratio due to application of a 
pressure of P = 0.75 GPa. 

the molecular model. We follow Ref. [13 in the use of the 
geometrical formulas t — {t2 + tn) /2 and t' — for 
obtaining t and t' of Eq. [TJ In order to corroborate these 
results, we use the alternative method of considering the 
ET dimers as sites of the tight binding model as shown 
in Fig.[T](d). Then the two dimers per unit cell for k-CN 
and K-SCN are responsible for the two antibonding ET 
bands at the Fermi level (or four dimers/four antibonding 
bands in the case of K-Br and k-CI) . 

The result of these extensive bandstructure calcula- 
tions and tight binding fits is summarized in Fig. ID The 
choice of basis set (LAPW or FPLO), functional (GGA or 
LDA) and number of sites included in the fit (molecules 
or dimers) results in a certain spread of resulting t, t' 
pairs which we interpret as a margin of error. The slopes 
of the lines connecting the orig in with the tight binding 
results for the GGA bands [31| show the results for the 
t' /t ratios at ambient pressure: t' /t = 0.83 ± 0.08 for n- 
CN, t'/t = 0.58 ± 0.05 for k-SCN, t' /t = 0.44 ± 0.05 for 
K-Cl, and t'/t = 0.42 ± 0.08 for K-Br. 

With the tight binding results, we can estimate the 
Hubbard [/ as C/ « 2ti following Ref. [13]. We obtain 
U/t ratios of U/t = 7.3 for k-CN, U/t = 6.0 for k-SCN, 
U/t = 5.5 for K-Cl and U/t = 5.1 for k-Bt. These val- 
ues provide a readjustment of the position of the four 
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K-type CT salts in the {U/t,t' /t) phase diagram, com- 
pared to the semiempirical calculations: For k-CN, we 
replace (8.2, 1.06) by (7.3,0.83), for k-SCN (6.8,0.84) by 
(6.0,0.58), for K-Cl (7.5,0.75) by (5.5,0.44) and for n- 
Br (7.2,0.68) by (5.1,0.42). At P = 0.75 GPa, we have 
(6.0,0.75) for K-CN and (5.7,0.58) for k-SCN. 

In Fig. |4] we also show the analysis of pressure ef- 
fects on the t, t' parameters. We represented with ar- 
rows the parameter changes from the ambient pressure 
results to the P = 0.75 GPa results in k-SCN and k- 
CN. For K-CN, pressure leads to a decrease in the t' /t 
ratio, while we observe almost constant t' /t for k-SCN. 
On the other hand we can establish that the parame- 
ter differences among the three studied cases at ambient 
pressure can be attributed to internal chemical pressure 
due to the different anion sizes. By comparing the ratios 
t' /t for K-CN, K-SCN and k-C1 at ambient pressure with 
the physical pressure t' /t behavior for k-CN and k-SCN 
we observe that chemical pressure and physical pressure 
have the same effect on the t' /t ratio: for k-CN, pressure 
significantly decreases the t' /t ratio. This effect is weaker 
for K-SCN. 

We now consider the new positioning of the four re- 
type CT salts in some of the calculated {U/t,t' /t) phase 
diagrams. In the PIRG phase diagram of Ref. IH, k-CI are 
shifted from the nonmagnetic insulating (NMI) phase to 
the antiferromagnetic insulating (AFI) phase which is an 
improvement as the superconducting phase is not con- 
sidered in this work. Meanwhile, k-CN and k-SCN are 
only repositioned within the NMI phase. In the CDMFT 
phase diagram of Ref. |3, K-Br is moved from the AF 
phase to the d-wave superconducting (SC) phase; k-C1 
is moved from the spin liquid (SL) phase to the border 
between AF and d-wave superconducting (SC) phase, k- 
CN is moved from the d-wave SC phase to the border 
between SL and SC phases, and k-SCN is moved within 
the SC phase. The effect of pressure means that k-CN is 
moved diagonally across the SC phase by the application 
of P = 0.75 GPa. In particular, the positions of K-Br, n- 
Cl and K-CN in the phase diagram are clearly improved. 
Our results exclude the possibility of viewing k-CN as a 
quasi-one-dimensional system [33| . 

In conclusion, we have presented the results of den- 
sity functional theory calculations for four K-type charge 
transfer salts. We obtain important shifts of the familiar 
t' /t ratios from semiempirical electronic structure calcu- 
lations [l3| towards significantly smaller values, i.e. to- 
wards lower frustration. Due to the smaller overall band 
width of the ET bands at the Fermi level, our estimate 
for the U /t values is also below the prevalent values. Our 
results call for a reexamination of the description of the 
K-type charge transfer salts, and in particular for k-CN, 
namely which is the nature of the spin liquid at ambi- 
ent pressure as well as the phase transition from Mott 
insulator to superconducting state under pressure. 

Note added: While finalizing this manuscript, Naka- 



mura et al. [3J| have posted a manuscript, where similar 
values oft'/t for k-CN and k-SCN as in our full-potential 
calculations are obtained from a pseudopotential method. 
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